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Display cases in supermarket systems often exhibit synchronization, in which the expansion valves 
in the display cases turn on and off at exactly the same time. The study of the influence of switching 
noise on synchronization in supermarket refrigeration systems is the subject matter of this work. For 
this purpose, we model it as a hybrid system, for which synchronization corresponds to a periodic 
trajectory. Subsequently, we investigate the influence of switching noise. We develop a statistical 
method for computing an intensity function, which measures how often the refrigeration system stays 
synchronized. By analyzing the intensity, we conclude that the increase in measurement uncertainty 
yields the decrease at the prevalence of synchronization. 


1 Introduction 

In a supermarket refrigeration, the interaction between temperature controllers leads to a synchronization 
of the display cases in which the expansion valves in the display cases turn on at the same time. This 
phenomenon causes high wear of compressors @. 

In this article, we apply the concepts in 0 to study a system in fl2l . More precisely, we combine the 
hybrid system model of a refrigeration system presented in llT2l with the proposed method for modeling 
switching noise presented in 0, and show that this formalism is able to characterize de-synchronizing 
behavior as a consequence of inaccuracy of temperature measurement. The affect of this inaccuracy has 
been utilized in a patent that proposes to adjust the cut-in and cut-out temperatures for the refrigeration 
entities to de-synchronize them lITOl . 

The definitions of synchronization have been formulated in 0, and numerous examples of synchro¬ 
nization have been analyzed in 0. In Ifl2l . the synchronization of a supermarket refrigeration system 
was studied as periodic trajectories of a hybrid system with a state space consisting of a disjoint union 
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of polyhedral sets and discrete transitions realized by reset maps defined on the facets of the polyhedral 
sets. 

However in real life applications, the switching instances are not deterministic, they are influenced by 
system and measurement noise. Therefore, we suggest a practical method rooted in stochastic analysis 
ItTTl . where noise is introduced into the system. For the stochastic analysis, we “thicken” the switch¬ 
ing surfaces to a neighborhood around them and formulate a probability measure of the switching in 
such a way that the longer the trajectory stays within the neighborhood, the higher is the probability 
of switching. As a result, this method provides the means of computing intensity plots. By analyzing 
them, we conclude that contrarily to the deterministic behavior, the system synchronizes intermittently. 
In addition, we infer that the prevalence of synchronization increases with decreasing uncertainties of 
the temperature measurements. In conclusion, the method of intensity plots provides a tangible test of 
whether and how frequently synchronization occurs in a refrigeration system. 

The article is organized as follows. In Section [2] and Section [3] we recall, from lfl2l . a simple model 
of a refrigeration system and demonstrate that this model coincide with the notion of a hybrid system, 
technical details are given in the appendix. Section [4] start by introducing the notion of switching un- 
certainity from 0, which then is used to incorporate inaccuracy of temperature measurement into the 
model. Subsequently a synchronization analysis of the refrigeration system is performed. 

2 Refrigeration system 

First, a brief description of a refrigeration cycle of a supermarket refrigeration systems with display cases 
and compressors connected in parallel, see figure Q]for a graphic layout. 

The compressors, which maintain the flow of refrigerant, compress refrigerant drained from the 
suction manifold. Subsequently, the refrigerant passes through the condenser and flows into the liquid 
manifold. Each display case is equipped with an expansion valve, through which the refrigerant flows 
into the evaporator in the display case. In the evaporator, the refrigerant absorbs heat from the foodstuffs. 
As a result, it changes its phase from liquid to gas. Finally, the vaporized refrigerant flows back into the 
suction manifold. 

In a typical supermarket refrigeration system, the temperature in each display case is controlled by 
a hysteresis controller that opens the expansion valve when the air temperature T (measured near to the 
foodstuffs) reaches a predefined upper temperature limit T u . The valve stays open until T decreases to 
the lower temperature limit T 1 . At this point, the controller closes the valve again. Practice reveals that 
if the display cases are similar, the hysteresis controllers have tendency to synchronize the display cases 
l6l . It means that the air temperatures 7} for i € {1,... ,1V}, where N is the number of display cases, tend 
to match as time progresses. 

In the sequel we discuss, for simplicity, a model of a refrigeration system that consists of only 
two identical display cases and a compressor. The dynamics of the air temperature 7} for display case 
i £ {1,2} and the suction pressure P for the system of two display cases are governed by the following 
system of equations, 


x= &(*) =A s x + B s , (1) 

where 
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Figure 1: Layout of a simple supermarket refrigeration system. 


with a,/3 constants whose specific values are provided in Appendix 16. 11 and 8 = (Si, 82) € 2 2 

with 2 = {0,1} and 5,- € {0,1} the switching parameter for display case 7; it indicates whether the 
expansion valve is closed (5, = 0) or open (t), = 1). The switching law is given by the hysteresis control: 


Si 


'l if 7] > 7f 
< 0 if 7 ; < T\ 

8i if T} < 7} < T‘‘. 


( 2 ) 


where T‘‘ and Tj are respectively the predefined upper and lower temperature limits for display case i. 
By convention, 5, = 0 for any initial condition 7j(fo) = 7j° e]7j ; , T-‘[. Such an initial condition is assumed 
throughout this paper; hence, (0 is well defined. 


3 Refrgeration system as a hybrid system 

Consider the following scenario. Let x(to) € ] T[ .7}'[x]T }.Tf [xl h and 8 = (0,0); thereby, both 
display cases are initially switched off. Suppose that at time t, the air temperature 7) of the 7th display 
case reaches the upper temperature limit T} 1 , then the 7th display case is switched on, and <5, = 1. This 
scenario indicates that the refrigeration system dT]) comprises four dynamical systems, each defined on a 
copy of the polyhedral set 

q = [t{,t?] x [r 2 ',r 2 "]xi +I (3) 

as illustrated in Fig.|2j A discrete transition between these four systems takes place whenever a trajectory 
reaches one of the following four facets of Q 


K = F\{Q) = S a [T’,Tl l ] x [Tin] x M+ and 
7f = F 2 “(G) = [T/,77] x 8 a [T‘,T?} x R+, 
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Figure 2: The state space of the refrigeration system consisting of two display cases is illustrated. Flere, 
the pressure axis is suppressed, and (Tj,T“) = (0,5). The direction of the vector field is indicated by 
the dark shaded triangles (see Proposition 1 in Ifl2l0 . 


with ae2 and 8 °[a,b] = {«}, 8 l [a,b\ = {b}. 

Moreover, the transitions between subsystems can be describe by eight reset maps Rj(8) : /r /,(, <5 ) x 
{5} x {/(/, <5)} defined by Rj(8)(x,8) = (. x,l(i, 8 )) with 


l{i, 8 ) = (h(i, 8 ),l 2 (i, 8 )) 


({81 + 1 , 82 ) if i= 1 
{( 81,82 + 1) if i = 2, 


where the results of the summation are computed modulo 2. Intuitively, the map I takes a polyhedral set 
enumerated by 8 to the future polyhedral set. The variable i indicates that the discrete transition takes 
place when the temperature 7} reaches its upper or lower boundary. The domain of a reset map, will be 
referred to as a switching surface. Specifically, a switching surface is a facet of one of the polyhedral 
sets making up the state space of the refrigeration system (QQ). 

The above construction allows us to identify the refrigeration system as a switched hybrid system. 
More precisely, let & = {Ps \Ps = Q x{5}, 8 G 2 2 } consist of the four polyhedral sets making up 
the state space, 5? = {^5 | 8,s( x ) = 4 gx + Bg, 8 G 2 2 } consist of the four dynamical systems, and 
= {Ri(8) | (i,S) £ {1,2} x 2 2 } consist of the eight reset maps. Then the triple (&,+?,&) constitute 
a hybrid system as defined in Appendix [672] Hence the behavior of the system is expressed by means of 
a (hybrid) trajectory which heuristically can be described simply as the union of trajectories generated 
by four local dynamical systems, see Appendix 16.3 1 for a precise describtion. 

With this notion at hand we say that a refrigeration system exhibit asymptotic synchronization if 
there exists a (T,l )-periodic trajectory which is asymptotically stable in A* [5] Definition 13.3], where 
X* denote the quotient space X/ ~ with X = Uge2 2 ^<5 ar *d ~C X xX the equivalence relation, O, 
generated by the reset maps in For a detailed explanation see Ifl2l . where it is also shown that the 
refrigeration system generates an asymptotically stable (T p , 2)-periodic trajectory, with T p ~ 261, lying 
on the diagonal of /} 0 , 0 ) and P, M . 

The (hybrid) refrigeration system with two display cases is illustrated in Fig. [3] Here, each element 
of & has been (orthogonally) projected onto the (7). 7?)-space. Hence, the polyhedral sets P§ are rep¬ 
resented by cubes. The three cubes 7} 0 ,t)> P(i,o)> P 1.11 have been vertically and/or horizontally reflected 
(compare with Fig. [2]). The stippled lines in the drawing indicate the reset maps in 
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Figure 3: The T\ TVstatc space of the refrigeration system consists of two display cases. The reset maps 
are indicated by the stippled lines (see Fig. |2] and its caption for further explanation). The pressure axis 
has been suppressed; thus, each P§ = Q x {5} is illustrated by a square. By abuse of notation, the facets 
of P§ are denoted by F[ x (instead of F[ x x {5})- 

4 Stochastic analysis of synchronisation 

So far, the discrete transition from a local system to another has been described as deterministic. In other 
words, it take place with probability one when a trajectory “hits” a switching surface. However, from a 
practical point of view, this causes a problem. For instance, any sampling step will result in the trajectory 
“hitting” the switching surface with probability zero. On account of model uncertainty, noise, and most 
of all the inaccuracy of temperature measurement, the switching surface should therefore be “thickened” 
by replacing each switching surface with an open neighborhood of it. 

In the sequel, we incorporate the “thickening of switching surfaces” into our model of the refrig¬ 
eration system. Specifically, we construct an open neighborhood around the switching surfaces within 
which a probability measure on each trajectory is proposed. Afterwards, we use this measure to describe 
the probability of a discrete transition, in such a way that the longer a trajectory stays within the neigh¬ 
borhood the higher becomes the probability of a transition. This makes it possible to devise a method for 
analyzing the typical behavior of the system, as well as a way of visualizing it. 

Let U be a stochastic variable uniformly distributed on [—e, e] with density function pin) = I[-e )£ ] («)/2e, 
where Ia denotes the indicator function of a set A. Consequently, the distribution of U can be described 
by the survivor function 

1 u < — £ 

/ U 

p{v)dv=< l — (u + e)/2e \u\ < e 

(O U>£, 

which is the probability ofU>u, P (U > u). The distribution of U can also be described by its conditional 
intensity (or hazard) function h{u) = p{u)/S{u) = I[_ e , e ] («)/(£ — u). Heuristically, a small variation of 
h(u ) is the probability of U being in a small region around u conditional on U not being smaller than u, 

P (u e du | u > u). The conditional intensity function, for short intensity, turns out to be a convenient 
starting point for defining stochastic transitions between the different local dynamical systems. 
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To this end, we derive an intensity function suitable for our puipose. Recall the notation introduced 
in Section I(x2l and let F(8) = f' 1 >1 ' 1 ' 5 ' 1 Uiv^' 5> denote the union of switching surfaces in the polyhedral 
set Pg. Furthermore, let F(8) e denote the £-neighborhood of F(8) in R 3 . 

For v e F(8) e , let be a switching surface within £ (Hausdorff) distance to x, and define 

x (i,S) =x ~ft(i,8)( x ), where is the orthogonal projection onto F Hence, when x (r g is non-zero, 
it is a normal vector to f ] i< ' 1 ' 5 ' 1 (more precisely, to the affine hull of F ^ (l8) ). It points into P§ when x g 
Pg—Ff'^’ S ^ and out of Pg when x € F(8) £ — Pg. Let »(,• §) denote a normal vector to F^ 1 ’ 8 ^ which points 
out of Pg. By means of the above quantities, we define the parameter U(- L g) (x) = |sign(/i ( , ^,v (; ^) 
where = x^g) (x) is regarded as a function of x. 

Now let y{t\k) denote a trajectory (see Appendix 16.31) of the refrigeration system, and assume that 
y(t\k) follows system 8, i.e., y(t‘,k) = E,g{y(t\k)). The intensity function h(ig) f° r switching from system 
8 to system l(i, 8) at y(t;k) can now be defined as 


h(i,5) (t(L&)) 


'o Uij(y(t;k)) < —£ 

< h(uij(y(t;k ))) \uij(y(f,k))\ < £ 

,°o Uij(y(t;k))> £. 


Hence, h^g^(y(p,k)) is the intensity of the signed orthogonal distance from y(t\k) to the switching sur¬ 
face F. !<lS> in the £-neighborhood of Pg. As a result, any trajectory whose intersection with the £- 
neighborhood of j s contained in a normal subspace to F- ,I ' l8> will switch according to the uniform 

distribution on the restriction of the trajectory to the £-neighborhood. 

The above construction successfully copes with switching for which a trajectory reaches one switch¬ 
ing region at a time. To deal with the situation where a point of switching is within distance £ to both 
switching surfaces, we will assume that switching to each of the systems will happen independently. 
With this assumption, we immediately conclude that if y(t;k) follows system 8, and is within distance 
£ from both F^~' 5 ^ and F^~' 8 ' then the intensity function hg for switching from system 8 to system 
l(i,8) at y(t;k) is given by hg(y(t\k)(t)) = h(\ g\{y{p,k)) + ftp, 5 ) (/(?;&))• Thus, if the discrete transition 
occurs at y(t',k), the switch to the system l(i,8) happens with probability h(^^{y{t\k))/h^{y{f,k)). 
Equivalently, we can allow switching to happen according to all of the intensities ,g)(y(t‘,k)) and 
h( 2 , 8 )(y(f,k)) independently of each other and disregard all the switchings except the first one which 
determines the switching. 

Based on the above theory and numerical simulations, we now conduct synchronization analysis of 
the refrigeration system when switching noise is present. To begin with, we describe how the theory is 
implemented in simulation. For this, we follow Algorithm 7.4III in {3[ p. 260]. 

When a trajectory enters the switching region F(8) e , an exponentially distributed variable with mean 
one is simulated. Thereafter, for each sample, the intensity is calculated, and subsequently the integral 
of the intensity is computed inductively. Afterwards, the exponentially distributed variable is compared 
with the number obtained by the integral computation [3] p. 258 (Lemma 7.4II)]. As a consequence, the 
switching occurs at the sample point where the integral exceed the exponentially distributed variable. 

In the above setup, the system trajectories become stochastic; hence, we can apply concepts such as 
mean values to describe the system behavior. For a trajectory (y, &*>), A C X and / = t',t"\ c [0, °°), 
let Zy(A,I) be the arc length from t' to t" of y inside A. As a consequence, y i-> Zy(AJ) defines a non¬ 
negative random variable, whose mean E[Z(A,7j] describes the curve intensity relative to (A,I). The 
curve intensity measures the typical behavior of trajectories by the mean length of trajectories on I inside 
A. 
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We illustrate the typical behavior of the refrigeration system with two display cases by simulation of 
the curve intensity. For this purpose, we use the following algorithm: (1) Fix a time interval I C [0 •> °°) 
and a sufficiently regular subset A C X. (2) Simulate n trajectories yi,...y„. (3) Divide A into small 
subsets Aj . (4) Calculate E[Z(A/,/)] ~ ^L/Zy^A,-,/) where Z-, ; /4,./) is approximated by the number of 
sampling points (on trajectories) falling in A,-. 



Figure 4: Each of the plots in column one and two are illustrated on various projections of the space 
obtained by identifying the four copies of Q in © used to define the state space of the refrigeration 
system. The first column represents intensity plots, £ = 0.1, I = [0,10 6 ], based on n = 10 trajectories at 
the point (0,0,15.23) on the asymptotically stable (T p , 2)-periodic trajectory found in Itl2l (top: 7) 73- 
plane, bottom: T\ P-plane). The second column represents one of the ten trajectories in a subinterval of 
I (top: T\ 73-plane, bottom: 7) P-planc) shown in the time interval [ 10 6 ,1.1 • 10 6 ]. The last two columns 
represent four of the ten simulations and show the sample-time evolution of suction pressure P. 

Intensity plots for I = [0,10 6 ], £ = 0.1 and n = 10 have been generated for initial values both within 
and outside the basin of attraction of the asymptotically stable (T pi 2)-periodic trajectory found in Ifl2l . 
The resulting intensity plots are similar to the one illustrated in the first column of Figure 0 Contrarily to 
the deterministic behavior, the system synchronizes occasionally over the time interval I, as seen in the 
last two columns of Figured where the pressure peaks at 15.23 correspond to the presence of the (T p .2)- 
periodic trajectory (synchronization). The behavior of the system in the remaining time is indicated in 
the first column of Figure 0 as the colored area except for the diagonal (top) and yellow/green circle 
through (0,15.23), and in the last two columns as pressure peaks at 10. We remark that the peaks at 
10 do not correspond to another limit cycle, as shown in the second column of Figure 0 (in fact this 
correspond to the interval [10 6 ,1.1 10 6 ] of the plot in the right lower corner). 

Further simulations reveal that the time spent in the (T p , 2) -periodic trajectory is highly dependent 
on £. For e = 0.1, relatively little time is spent in the (T p , 2)-periodic trajectory, whereas for £ = 0.01, 
the behavior closely matches that of the deterministic case. More precisely, for £ = 0.01 and initial 
values in the basin of attraction, the system rapidly converges to the (T p , 2) -periodic trajectory and stays 
close to it with a high probability, while for initial values outside the basin of attraction, the system 
has a small probability of converging to the (T p , 2)-periodic trajectory within the time interval 1. As 
a consequence, the relative frequency of appearance of the (T p , 2)- trajectory can be described as a 
monotonically decreasing function /(e) € [0,1], where /(0) = 1 corresponds to the deterministic case. 

In summary, uncertainties of the temperature measurements influence synchronization and can be 
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used to design de-synchronization scheme. We conclude that adding noise as above to temperature 
measurements de-synchronizes the refrigeration system: An increase in measurement uncertainty yields 
a decrease of the prevalence of synchronization. 


5 Conclusion 

In this paper, we have investigate the influence of switching noise on the synchronization phenomenon 
in supermarket refrigeration systems. We have developed a numerical method for computing intensity 
plots. By analyzing them, we have concluded that the time the refrigeration system stays synchronized is 
dependent on uncertainties of the temperature measurements. A greater measurement uncertainty yields 
a smaller accumulated time in which the refrigeration system is in synchronization. 


6 Appendix 

6.1 Model of refrigeration system 

The mathematical model presented here is a summary of the model developed in fTTll . For the ith display 
case, dynamics of the air temperature T a - ir i can be formulated as 


dT m 


Qgoods—air,i T Qload,i &iQe,i 


,max,i 


dt 


1 + 


UA, 


goods—air, i 


UA a if— wa ll 

Qgoods 


with 


’. j ^wall,i^p,wall,i 
air.i Qload.i 


Twallj — Tair,i j, * 7 

U A-air—wall,i 

Qgoods—air, i = UA. goods—air,i{TgO,i T a ir,i ) ? 

Qe,max,i — UA wa ii— re f ,max,i ( .i O-jPsUC k'f ) ; 


(4a) 

(4b) 

(4c) 

(4d) 


where the process parameters arc specified in Table [Q and 5, G {0,1} is the switch parameter for the ith 
display case. When 5,- = 0, the ith expansion valve is switched off, whereas when 5,- = 1 it is switched 
on. The suction manifold dynamics is governed by the differential equation 

dPsuc SiTh 0 ,i T dl r const V comp iP’P Pme + b p ) 

dt Vsuc ' ^ PsucO 

where N is the number of display cases, and mo,; = /ho for i £ {1,... ,N}. 

We denote 7} = 7)„y ; and P = P suc and write the dynamics of the air temperature and suction pressure 
in the concise form (with the process constants in (l4al) collected in a.b,c,d,e. a,[5 and then replaced by 
their numerical values) 


Tj = —aTj + b — Si ( cTj — dP — e) (6a) 

= -0.00197)- + 0.0244 - 5, (-0.00127)- + 0.0506P + 0.1065), 

P = -aP + f5 + 8 l + S 2 = —0.056 P + 0.0038 + 5i + S 2 . 


(6b) 
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Table 1: Parameters for a simplified supermarket refrigeration system 

6.2 Hybrid systems 

A detailed study of the hybrid system presented below can be found in 0. We write F A P if F is a face 
of the polyhedral set P. A map / : P —> P' is polyhedral if 1) it is a continuous injection, and 2) for any 
F yP there is F' -< P' with dim(F) = dim(F') such that f(F ) C F'. 

Definition 1 (Hybrid System) For finite index sets J and D, a hybrid system (of dimension n) is a triple 
,3/ = (SPrhdf ifi&j), where 

1. & = {P S C M" | P$ a polyhedral set. d i m (P§) = n. 8 € D } is a family of polyhedral sets. 

2. SP = K" | P§ G ffd ■ 8 GO} is a family of smooth vector fields. 

3. .% = {Rj :F -+F' \F -<P G F' -<P' G dim(F) = dim(F') =n—l, ;G/} is a family of 
polyhedral maps, called reset maps. 

After identifying D with a finite subset of M, we can rewrite the hybrid system {fido,//' ,Mj) as Fig. 

e F(x,q) = {($ 9 (x),0)| x G P q } for (x,q) G C, and 
(x,q) + G {Rk(x,q) \ x G G(x,q ) = dom(/4) -< P q } for (x,q) G D , 

where C = U Sgd( p 8 x {5}) and D = U {(8j)eDxj | dom (Rj)^p 5 }( dom (Rj) x {5}). This is precisely the 
hybrid system in | f 4j|. 

6.3 Trajectories of a refrigeration system 

We bring in a concept of a (hybrid) time domain 0|. Let k G N U {°°}; a subset C R + x Z + will be 
called a time domain if there exists an increasing sequence {hhe { 0 ,,,,,/c} i n ® - U {°°} such that 

= U P x « 

where T { = [»/_i, t t \ if i G { 1,..., k - 1 }, and T k = j ^ < “ 

Note that T, = for all i if k = °o. We say that the time domain is infinite if k = °° or 4 = °o. 

The sequence {4 } ze {o,.._ ,*} corresponding to a time domain will be called a switching sequence. 

Definition 2 (Trajectory) A trajectory of the hybrid system (/jd n .y' fff.j) is a pair ( ■ 7) where k G 
N U {°°} is fixed, and 
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• 3F k C M + x Z+ is a time domain with corresponding switching sequence {fi}, G {o,. 

• 7 : l? k —>X = x {^i / ' v continuous (X has the disjoint union topology) and satisfies: 

1. For each i £ {1 ,... ,k— 1}, there exist S / 3' £ D swcfi that y(tf, i) £ bdjPg), and y(t,U+ I) £ 
bd(P§i). 

2. For each i £ { I..... A:}, there exists 8 £D such that the Cauchy problem ijit-, 1 ) = y{t-,i) = 

<^5(7(4;/)), y(ti-\,i) = £ Pg, lias a solution on Tj C LP k . 

3. For each i £ {1, 1}, there exists j £ J such that Rj(y(tf,i)) = 7(4,;/+1). 

A trajectory atx is a trajectory (,%. y) with y(to', 1) = x. By abuse of notation ywill sometimes be referred 
to as a trajectory. 

The next definition formalizes the notion of a periodic trajectory, which will be used in defining 
synchronization of the refrigeration system. 

Definition 3 ((T,/)-periodic trajectory) Let (T,l) £ M+ x Z + . A trajectory (■%. y) is (TJ)-periodic 
(or just periodic) if(l) LP k is an infinite time domain, and (2) for any i £ {1,... ,k} and t £ p(£L k ), where 
p : FL k —> [?o, °°[ is the projection p(t , i ) = t, we have y(t + T:i + l) = 7(4; i ). 

In particular, if (■%, y) is a (P,/)-periodic trajectory, and T is nonzero then p : ■% [to, 00 [is surjec¬ 
tive. 
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